Dynamic critical behavior of model A in films: Zero-mode boundary conditions and expansion near 

four dimensions 
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The critical dynamics of relaxational stochastic models with nonconserved 7i-component order parameter 
(p and no coupling to other slow variables ("model A") is investigated in film geometries for the cases of 
periodic and free boundary conditions. The Hamiltonian Tt governing the stationary equilibrium distribution 
is taken to be 0(n) symmetric and to involve, in the case of free boundary conditions, the boundary terms 



/2 associated with the two confining surface planes *8j, j = 1, 2, at z = and z = L. Both 



enhancement variables Cj are presumed to be subcritical or critical, so that no long-range surface order can 
occur above the bulk critical temperature r c .oo. A field-theoretic renormalization-group study of the dynamic 
critical behavior at d — 4 — e bulk dimensions is presented, with special attention paid to the cases where 
the classical theories involve zero modes at T Ci00 . This applies when either both Cj take the critical value c sp 
associated with the special surface transition or else periodic boundary conditions are imposed. Owing to the 
zero modes, the e expansion becomes ill-defined at T Cj00 . Analogously to the static case, the field theory can 
be reorganized to obtain a well-defined small-e expansion involving half-integer powers of e, modulated by 
powers of In e. This is achieved through the construction of an effective (d — l)-dimensional action for the zero- 
mode component of the order parameter by integrating out its orthogonal component via renormalization-group 
improved perturbation theory. Explicit results for the scaling functions of temperature-dependent finite-size 
susceptibilities at temperatures T > T Ci00 and of layer and surface susceptibilities at the bulk critical point are 
given to orders e and e 3 ^ 2 , respectively. They show that L dependent shifts of the multicritical special point occur 
along the temperature and enhancement axes. For the case of periodic boundary conditions, the consistency of 
the expansions to 0(e 3 ^ 2 ) with exact large-n results is shown. We also discuss briefly the effects of weak 
anisotropy, relating theories whose Hamiltonian involves a generalized square gradient term B kl dk<fr ■ di<p to 
those with a conventional (V0) 2 term. 



I. INTRODUCTION 

The renormalization group (RG) approach has played an 
important role in the modern theory of critical phenomena. 1 "* 
For one thing, it provides an appropriate mathematical frame- 
work for the formulation of the theory. Second, it has led 
to the development of powerful calculational tools for quan- 
titatively accurate investigations. Its most impressive and nu- 
merous successes have been achieved in the study of static 
bulk critical phenomena. However, appropriate extensions for 
studies of dynamic bulk critical phenomena, 5,6 boundary crit- 
ical phenomena, 7 8 and finite-size effects 9-22 were developed, 
which have proven their utility and power. Anyone of the fea- 
tures, "dynamics", "boundaries", and "finite size," involves 
fundamental new issues and adds to the technical complex- 
ity of analytic RG studies. It is therefore not surprising that 
work on problems involving combinations of several of these 
features has remained rather scarce. 

In this paper we shall be concerned with the dynamics of 
systems in slabs R d_1 x [0, L] of finite thickness L near then- 
bulk (L = oo) critical point. Both the cases of periodic and 
free boundary conditions will be considered. Hence we shall 
have to deal with all three of the above-mentioned features. 

Despite the availability of some experimental results, there 
exist only few previous studies of dynamic critical behavior of 
systems in film geometry (see Ref. 23 and its references). In 
two earlier papers, Calvo and Ferrell 24 25 investigated the dy- 
namics of binary liquid mixtures confined between two par- 



allel plates using the mode-mode coupling approach. Sub- 
sequently, the dynamics of bounded one-component fluids 
near the liquid-gas critical point (model H in the terminol- 
ogy of Ref. 5) and of confined liquid 4 He at the super- 
fluid transition (model E) were analyzed by this method. 26 
There exist also a number of papers dealing with finite-size 
effects on dynamic critical behavior and dynamic surface 
critical behavior. 18 ~ 20 ' 22,27 ^ 3 In other works, finite-size sys- 
tems with long-range interactions or quenched disorder were 
investigated. 44 ^ 6 

Recently, Gambassi and Dietrich 47 presented a fairly de- 
tailed study of the familiar model A (Refs. 5 and 6) in film 
geometry within the framework of the classical (zero-loop, 
van Hove) approximation, augmented by RG-improved per- 
turbation theory. They focused on the situation where the sur- 
face interactions on both confining surfaces are subcritically 
enhanced. This corresponds to the case in which ordinary sur- 
face transitions 7,48 occur in the semi-infinite systems bounded 
by either one of the two surface planes. Knowing that Dirich- 
let boundary conditions apply under these conditions at both 
boundary planes on sufficiently large length scales, they re- 
stricted their analysis by choosing such boundary conditions 
from the outset. 

In the present paper we shall also be concerned with model 
A in film geometry. Our analysis complements and goes be- 
yond that of Ref. 47 in several ways. First, we shall not 
limit ourselves to the classical approximation but present a RG 
analysis in d = 4 — e bulk dimensions, going to one-loop order 
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in our explicit calculations (and partly beyond to determine 
contributions of order e 3 ' 2 ). Second, we shall give up the 
restriction to Dirichlet boundary conditions on both confin- 
ing plates. Aside from periodic boundary conditions, we shall 
consider, in the general part of our analysis, the generic case of 
symmetry preserving Robin boundary conditions correspond- 
ing to distinct enhancements of the surface interactions on the 
two boundary planes. This includes the case of special-special 
(sp-sp) boundary conditions for which the surface interactions 
on both boundary planes are critically enhanced. 49-51 



Periodic and sp-sp boundary conditions share the feature 
that Landau theory involves a zero mode at bulk criticality. 
It has recently become clear that this causes a breakdown of 
the e expansion at T c>00 . 50,5 The small-e expansions of the 
associated universal amplitudes of the critical Casimir forces 
were found to involve, besides integer powers of e, also frac- 
tional powers e k < 2 with k > 3 (modulo powers of In e). This 
breakdown of the e expansion is similar to the one reported in 
Ref. 52 for the 3 — d expansion of bosonic quantum systems. 



II. THE MODEL AND BACKGROUND 

A. Definition of model A in film geometry 

We begin by defining an appropriate extension of model A 
for the film geometry. To this end, we consider a film occupy- 



ing the region 23 = 



pd-l 



[0, L] of d-dimensional space 



R d . We wish to study the critical dynamics of such films 
involving an 77-component order-parameter field <p(x, t) = 
(<j> a (x, t),a = 1, . . . , n). We write position vectors as x = 
(y, z), where y G W 1 ^ 1 and z G [0, L] are the coordinates 
alongside and across the film, respectively. We choose peri- 
odic boundary conditions along the d — 1 principal y direc- 
tions. Depending on whether we are concerned with peri- 
odic or free boundary conditions in the z -direction, the slab 
23 has no boundary, <923 = 0, or consists of the two (d — 1)- 
dimensional confining hyperplanes 23 1 at z = and 232 at 
z = L. In the latter case, we orient the boundary such that the 
normal n on 523 = 23 = 23i U 232 points in the interior of 23. 

We are interested in the dynamics of systems that relax to a 
stationary equilibrium state described by the Hamiltonian (in 
units of ksT) 



The primary aim of this paper is to show that a similar 
breakdown of the e expansion is encountered in the study of 
dynamic critical behavior of model A in film geometry for pe- 
riodic and sp-sp boundary conditions. We shall demonstrate 
this explicitly by determining the contributions of order e 3 / 2 
of the dynamic finite-size susceptibility \l for both bound- 
ary conditions, the layer susceptibility xiz {L) for periodic 
boundary conditions, and their surface analogs Xii P sp ' (L) and 
ail (-k) f° r S P- S P boundary conditions at T c oo . 



The remainder of this paper is organized as follows. In 
Sec. II, we first define an appropriate extension of model A 
to the film geometry. We then recall the Lagrangian formu- 
lation of the corresponding Langevin equation 27 53 55 along 
with some necessary background such as the fluctuation- 
dissipation theorem and the boundary conditions of the order- 
parameter field <fr and the associated response field </>. 36,37 In 
Sec. Ill, we set up perturbation theory, explain the renormal- 
ization of the theory for d = 4— e, give the RG equations of the 
multi-point correlation and response functions, and describe 
their solutions. Section IV begins with a discussion of the ba- 
sis of RG-improved perturbation theory. Next, we show that 
the e expansion breaks down at T Cj00 for periodic and sp-sp 
boundary conditions and elucidate the origin of the problem. 
To obtain well-defined small-e expansions, we then construct 
an effective (d — l)-dimensional dynamic field theory for the 
zero-mode components of <f> and <f>. In Sec. V we present 
one-loop results for various scaling functions for T > T C)00 . 
Section VI contains a brief summary and concluding remarks. 
In addition, we briefly embark on the issue of universality vi- 
olations due to weak anisotropy and other sources brought up 
recently. 56 - 57 Finally, there are three appendixes in which tech- 
nical details are described. 



H\4>] = / d d x 



(2.1) 



Here the contributions localized on the boundary planes 23j, 
given in the second line of Eq. (2.1), are only present for 
free boundary conditions (p = /) but absent for periodic 
boundary conditions (p = per). The absence of bound- 
ary terms linear in <f> reflects our assumption that the bound- 
aries do not break the <j> — > — <fi of the Hamiltonian. That 
no quadratic anisotropies have been taken into account in 
the boundary terms meets our stronger requirement that the 
boundaries do not break the presumed 0(n) symmetry. Sur- 
face spin anisotropies, which would require separate enhance- 
ment variables c J Q for the boundary terms oc 2 of different 
components a, 58,59 will not be considered here. 

For the sake of simplicity, we shall furthermore assume 
that the values of both surface variables Cj are such that no 
long-range surface order can occur at 23j above the bulk crit- 
ical temperature T Ci00 . Recall that in a semi-infinite system 
bounded by 23 j, a transition to a bulk-disordered, surface- 
ordered phase takes place at a temperature T c s > T Ct00 when 
Cj drops below the threshold value c sp associated with the so- 
called special transition (provided the dimension d is suffi- 
ciently large that the d — 1 dimensional surface can support 
long-range order). Thus, our assumption translates into the 
conditions 



Scj = Cj - Csp > , .7 = 1,2 



(2.2) 



their physical meaning is that the surface pair interactions are 
subcritically (Scj > 0) or critically (S£j = 0) but not super- 
critically (Scj < 0) enhanced. 
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We are now ready to define an appropriate extension of 
model A to the film geometry considered that is compatible 
with our assumptions. Straightforward considerations analo- 
gous to those made in Refs. 36 and 37 for semi-infinite sys- 
tems lead us to consider the Langevin equations 



<p a (x,t) 



\ 6H t * 

0(Pa 



(a(x,t) 



(2.3a) 



which are meant in the sense of Ito. 60 Here A is the bare On- 
sager coefficient and £ is a Gaussian random force of mean 
zero, 

(Ca(a:,*)> = 0, (2.3b) 

and variance 

(Ca(a5,t)Ca'(aj',*')> =2\5 aa <5(x-x')5(t-t') . (2.3c) 

B. Lagrangian formulation of the theory 

In our subsequent analysis of this model it will be conve- 
nient to use its equivalent Lagrangian formulation. 27 - 36,37,53 ~ 55 
This involves the action 



A( V-V+f + -^)0 



0- A0 



S p.f ° x Y] / 



(2.4) 



where <fi is an auxiliary field, the so-called response field. The 
gradient operators V and V act as indicated to the left and 
right, respectively. 36,37 ' 61 Note that we have dropped a con- 
tribution oc 9{t = 0) [where 6(t) is the Heaviside function] 
produced by the Jacobian det(5^ a /64>p), choosing a prepoint 
discretization in time. 

We fix the initial condition for the solutions to Eq. (2.3) 
in the infinite past, taking the limits t{ — > — oo and 
tf — > oo, and suppress these integration limits henceforth. 
Multipoint correlation functions of the fields (f> a (x,t) and 
(f>a(x,t) can then be calculated with the functional weight 
cxp(— J\4>, 4>})V[<p, <p], where the measure V[<p, <p] is pro- 
portional to Y\ „ * d{(j> a /2m) d<p a and normalized such that 



T>[(t>, 4>] e~ J[ ^ ] = 1 



(2.5) 



The field <fr a (x,t) describes how averages (O [<p]) of ob- 
servables 0[4>] obtained from the solutions to the Langevin 
equation (2.3) upon averaging over noise histories respond 
to perturbations that change its right-hand side by a function 
J a (x, t). In the case of model A, the addition of the time- 
dependent magnetic-field terms 



/ d d xh(x,t) ■ 4>(x,t) 
J<n 

-8 pJ f cf^y/^OM) i) (2.6) 



to Hamiltonian (2.1) would yield such a perturbation with 
J a (x,t) = Xh a (x,t) force ^ 53 and corresponding bound- 
ary terms oc . To explain the consequences, let us introduce 
the generating functional 

g[J,K;J,K) = In^expij dt J (J-0 + J-0) 

(K-j> + K-<t>) |V (2.7) 

where K and K are source functions localized on the bound- 
ary ( B. They serve to generate the boundary operators (fy^ — 
4>(x<B,t) and </> s = <fi(x<B,t) with x<s e *B by functional 
differentiation. In the case of periodic boundary conditions, 
they are not needed, and we write Q[J;J] for the analog of 
functional (2.7). 

To specify a boundary point x<%, we must say on which 
boundary plane Q3j it is located and give its lateral coordinate 

r. Denoting the restriction of <f> a f to 93^ by <f> a 3 , we can write 

the boundary operators as (f>a 3 (r, t) and 0„ 3 (r,t). Whenever 
we do not wish to specify on which surface plane these bound- 
ary operators are localized, we continue writing </>® and <p^ . 
Analogous conventions will be used for the boundary sources 
Ka 3 and Kq 3 and the boundary magnetic fields A„ 3 . 
Functional (2.7) generates the cumulants 



N 



M 



6 >B 



M N 

\i=l k=l 1=1 m=l I 

EE W (N,M;N,M) (2 g) 

whose N + N + M + M position vectors, time arguments, 
and tensorial indices we have suppressed. From the corre- 
spondences 



XS 



SJ a (x,t) Sh a (x,t) 



^\4> a {x,t) (2.9) 



and 



6K**(r,t) 5h* 3 (r,t) 

it is clear that the functions wW^iv.Af) defined in Eq ( 2 .8) 
are the usual connected correlation and response functions. 



<-A^(r,t) (2.10) 



C. Fluctuation-dissipation theorem and mesoscopic boundary 
conditions 

Several other remarks are in order here, which concern the 
fluctuation-dissipation theorem and the boundary conditions 
on the mesoscopic scale where our continuum approximation 
applies. 

First, the fluctuation-dissipation theorem 5 



= -e(t-t')d t (Mx,t)^(x',t')) c 



(2.11) 
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holds. Second, the boundary contributions to the classical 
equations of motions yield the boundary conditions 



d n <f>(x, t) = c.j <t>{x, t) , x e 58j , 

d n 4>{x, t) = Cj 4>(x, t) , X e . 



(2.12) 



These hold beyond the classical approximation inside of aver- 
ages (up to anomalies at coinciding points). 36-38 On the level 
of the classical (zero-loop) approximation, they ensure that 
the matrix integral kernel 



2A 



9 t + A(-A + f) 



® (<W) 



<T(2) = 

' -d t + X(-A + r) 

(2.13) 

of the quadratic part of J[4>,4>] is self-adjoint. Third, 
the boundary conditions (2.12) imply that the fluctuation- 
dissipation theorem (2.1 1) remains valid when either x or x' 
approaches a surface point. 

In the case of periodic boundary conditions, we have 

4>(x + Le z ,t) = <p(x,t) , 

c/>(x + Le z ,t) = <j>(x,t), (2.14) 

instead of Eq. (2.12). 

III. PERTURBATION THEORY AND RG 

A. Free propagators 

To set up perturbation theory, we employ dimensional reg- 
ularization and focus on the disordered phase. The free re- 
sponse and correlation propagators, Rl and Cl, then follow 
from the inverse of the matrix kernel (2.13). We have 



J 



(2) 



( o 4 

\Rl Cl 



(3.1) 



where Rl is the solution to 



[d t + A(-A + f )] R L (x, t; x', if) = 5{x - x') 5(t - if) , 

(3.2) 

while Cl is proportional to the convolution Rl * R L : 

C L {x, t; x', t') = 2A (R L * R L )(x, t; x' , if) 

= 2A / dt I d d x R L (x,t;x,t) R L (x,t;x',t') .(3.3) 
J -co J<a 

The yi-Fourier transforms of these quantities (for which we 
use the notational conventions summarized in Appendix A) 
can be expressed as 



R L (p; z, z<- uj)=y fnWfW) 

-iu; + \(T + p2 + kl) 



(3.4) 



and 



C l (p;z,z';uj) = ^ 



2\f m {z)f* m {z') 



-iu + X(i+p 2 + kl l )\ 



(3.5) 



in terms of a complete set of orthonormal eigenfunctions 
fm{z) of the operator — d\, where fm(z') is the complex con- 
jugate of f m (z'). These functions are properly normalized 
solutions to 



&z fm(z) — k m f m {z) 



subject to the boundary conditions 



.4(0) = c 1 f m (o) , 

-f' m {L) = c 2 f m (L) • 



(3.6) 



(3.7) 



For non-negative values of c\ and 62, the spectrum {fc^} 
is discrete with kf n > 0. The eigenfunctions are phase- 
shifted cosine functions f m (z) = A m cos(k m z + $ m ), whose 
phase shift d m follows from the first of the boundary con- 
ditions (3.7). The eigenvalues are solutions to the transcen- 
dental equation implied by the boundary condition at z = L. 
For general values of Cj > 0, the eigenvalues kf n depend 
on both ci and 62, as do the normalization factors A m and 
the phase shifts d m (via fc m ). 62-64 For the special values 
(ci,C2) = (00,00), (0,0), and (00, 0) corresponding to the 
combinations D-D, N-N, and D-N of Dirichlet (D) and Neu- 
mann (N) boundary conditions on the two planes, the eigen- 
values kf n and eigenfunctions f m can be found in Appendix 
A of Ref. 49 and Appendix A of Ref. 78. 

However, the response propagator Rl(p', z, z';u>) can 
also be determined by solving the analog of Eq. (3.2) in 
the pzuj representation using familiar methods for Sturm- 
Liouville differential equations. 65 We give the result of such 
a calculation for general non-negative values of 61, C2, and 
f in Eq. (B5) of Appendix B. In the special cases t\ = 
c>2 = and b\ = 02 = 00, the result is equivalent to the 
representations 7 ' 10 ' 47 

( N " N ) 

R L D - D (p;zi,z 2 ;u}) = [ R °o(p, zi - z 2 - m2L, u) 

m— — 00 

± Rcoip, Z1 + Z2- m2L, cj)] (3.8) 



of the corresponding Neumann and Dirichlet propagators 
i?i N " N) and 4 D - D) as a sum of image contributions involving 
the bulk propagator (see, e.g., Refs. 7, 10, and 51) 



R 00 (p,z V2 ,uj) = T X- e -*«l»«l , 



2\k u 



where 212 = 21 — 2:2 and 



p 2 + t — iuj/X . 



(3.10) 



In the case of periodic boundary conditions, one has 

00 

R L pei \p;z 1 ,z 2 ;uj) = ^ i?oo(p, Z12 - mL, u) . (3.11) 



m— — 00 



The corresponding finite-size correlation propagators 
Cj j N N ' > , C L D D \ and Cjf 6 can be expressed in terms of the 
free bulk correlation propagator Coo(p, zi 2 ,u>) in a manner 
completely analogous to Eqs. (3.8) and (3.11). 
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B. Reparametrizations 

As is well known and explained elsewhere, 7 ' 10,20 - 27 - 36 ' 37 the 
singularities of -Roo(sEi2> ^12) ar, d Cao( x i2> £12) at coinciding 
points (#12, £12) = (0,0) produce ultraviolet (uv) singulari- 
ties in Feynman integrals of the multipoint cumulant and re- 
sponse functions (2.8). For dimensions d < 4 and periodic 
boundary conditions, the uv singularities of these functions 
can be absorbed via standard "bulk" reparametrizations of the 
form 



= 


Z <f, <PR, 


(3.12a) 


4> = 




(3.12b) 


A = 




(3.12c) 




\1 2 Z T T, 


(3.1 2d) 


uN d = 


H e Z u u. 


(3.12e) 



Sf = 



Here fi is an arbitrary momentum scale and f c denotes the crit- 
ical value of f of the d-dimensional bulk theory. Following 
Ref . 5 1 , we choose the factor that is absorbed in the renormal- 
ized coupling constant u as 



2T(3-d/2) 



2)(47r) d / 2 

1 - IE 



16vr 2 



1 



ln(47r) 



0(e 2 ) 



(3.12f) 



where 7^ = — r'(l) is Euler's constant. If we employ dimen- 
sional regularization and fix the renormalization factors Z g , 
g = <p, <fr, u, t, by minimal subtraction of poles in e, this con- 
vention ensures that the two-loop results for these functions 
given in equations (3.42a)-(3.42c) and (4.54) of Ref. 7 apply. 

For free boundary conditions, additional primitive uv sin- 
gularities with support on Q5i and S2 occur. These can be ab- 
sorbed through the additional ("surface") reparametrizations 



c j ~ c sp — (J-Z c Cj , 



(3.13) 



known from the semi-infinite case, where [0 ],r and [<j> }r 
are renormalized boundary operators. Explicit two-loop ex- 
pressions for the renormalization factors Z\ and Z c may be 
found in Eqs. (3.66a) and (3.66b) of Ref. 7 or in Refs. 66 and 
67. 



C. RG equations and scaling 

Upon introducing the renormalized functions 



W-. 



(N,M;N,M) 



R 



-(AM- AO/2 



(^Zi)"(* +M i/ a W&MW , (3.14) 



we can exploit the invariance of the bare functions 
^y(N,M;N,M) un( j er changes ji — > \il in a standard fashion 
to obtain the RG equations 



2? M + (#+#)^ + (M+M) 



w 



(N,M-N,M) 



R 



Here 



g—u,r,X.ci jC2 



. 
(3.15) 

(3.16) 



The beta and exponent functions (3 g and rj g are given by 

[3 g = fid^og = -{d g +T} g )g 
-(e + r) u )u, g = u, 
-(2 + t} t )t, g = r, 
(2 - ry A )A , g = X, 
-(1 + r} c )cj , g = Cj, j = 1,2. 



(3.17) 



and 



rjg = fid^log , g = u,r, <f>, <f>, A, 1, c-y, c 2 , 



(3.18) 



respectively, where d^o means a derivative at fixed parame- 
ters of the bare theory. Explicit results for (3 U to order u 3 and 
for n<j>(u), rj T {u), r}\{u), and r] c (u) = rj Cj (u) to order u 2 may 
be looked up in Eqs. (3.75a), (3.75b), (3.76a), and (3.76b) of 
Ref. 7, respectively. The function T]\, which can be written 
as rj\ = (77^ — rji)/2, is given to order u 2 in the first line of 
Eq. (III.8) of Ref. 27. 

It should be obvious how the above RG equations carry over 
to the case of periodic boundary conditions: the RG equation 

for W R N ' N) agrees with Eq. (3.15) if we set M = M = 0, 
except that the contributions to 2? M involving d Cj drop out be- 
cause these variables do not occur in the theory when periodic 
boundary conditions are chosen. 

Using characteristics the RG equations (3.15) can be ex- 
ploited in a familiar fashion to derive the asymptotic scaling 



behavior of the functions W 
to the flow equations 

dg{i) 



(N,M;N,M] 
R 



. Let g(£) be solutions 



£■ 



d( 



P g [u{£),g{t)}, g = u,T,X, Cl ,c 2 , (3.19) 



satisfying the initial conditions g(l) = g. Then we have 

dv 



ln£ = 



u Pu{v) 



(3.20) 



and 



g{t) = E g [u(£),u]£- ( - d ° +r >^ g, g = T,X,c 1 ,c a , (3.21) 



where E g (it, u) are the trajectory integrals 



E g (u, u) = exp 



dv 



(3.22) 
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and the asterisk indicates values rj* = at the nontrivial 

root u* = 0(e) of (3 u {u). 

To illustrate the consequences of the RG equa- 
tions (3.15), let us consider the multipoint cumu- 

lants in the cci-presentation. Writing w^< M < N > M ) _ 
Wr(x, t;u, t, Ci,c%,L, A, fi), where x and t represent all 
position and time variables, respectively, we note that the 
dimension of this function is \ N+M t with 

d w = (N + M) + (N + M) . (3.23) 

Using this in conjunction with the solution to the RG equa- 
tion (3.15), one arrives at 

W R {x, t; u, r, ci, c 2 , L, A, fx) = X N+M £ Aw 

xE w (£) W R {[iix, X/J.H; u, f , c 1 ,c 2 ,fj£L, 1,1). (3.24) 

Here is the scaling dimension of Wr. We use the notation 
A[0] for the scaling dimension of an operator 0(x, t). Then 
we have 

A[0] ee {d-2 + r$)/2 = p/v, 
A[$\ ee 2-^ + A[0] = j + /3/i/, 
A[^} ee (d-2 + r,; + r&)/2 = 0?/v, 
A[0 S ] ee t + A[<f>*]=t + [3?/v, (3.25) 

where we have introduced the standard static bulk critical 
indices j3 and v, the dynamic bulk critical exponent 3, and 
the surface critical exponent /3^ p of the special transition. 7 In 
terms of these quantities, the scaling dimension reads 



A w = NA[4>] + NA[<f>] + MA[ 



MA\ 



(3.26) 



Further, the prefactor Ew (£) denotes the trajectory integral 



Ew{£) = cxp 



dv 



(3.27) 



with 



. . N + N M + M . ~ 

Vw{u) = — - — v<f> H 2 W + fh) — \N + M) rj\ . 

(3.28) 

It is easily checked that Eq. (3.24) yields scaling forms in 
accordance with the phenomenological theory of finite-size 
scaling. 14 ' 15 Assuming that the initial coupling constant u is 
nonzero, we replace the running coupling constant u and the 
scale-dependent amplitude Ew{£) by their respective long- 
scale limits u(0) = u* and E^ = £w(0) and substitute the 
running variables f , A, and £j by their limiting forms 

f E*(u) l~ x l v T, 
g(£) w E* g (u)r^ + ^g = I E* x (u)P A, 

(3.29) 



where E*(u) = E g (u*,u) are nonuniversal amplitudes, and 
we have introduced the surface crossover exponent $ = + 
77*). 7 If we now fix the scale parameter £ = i T for given t^O 
by f(£ T ) = sign(r), then 



(3.30) 



agrees with the second-moment bulk correlation length 

£± « d 0) it/t C)0O - ir , 



(3.31) 



up to a nonuniversal scale factor. We thus obtain the scaling 
forms 

W R {x, t; u, t, ci, c 2 , L, A, /x) « A* +a7 / w 

x (M^)- Aw -B^H ± (x/e±,t;c 1 ,c 2j L/^±) (3.3,2) 

in which t represents the set of all scaled time variables 

t = ^V*0^i)" , ) (3-33) 

while ci and c 2 denote the scaled surface variables 

c, = £ c * Cj )* . (3.34) 

The scaling functions S± are given by the restrictions of 
Wr(. . . ; u, f, . . . ) to the hyperplanes with u = u* and 
f = ±1. 

Analogous results apply for periodic boundary conditions, 
except that the surface scaling variables ci and c 2 are missing. 



IV. RG-IMPROVED PERTURBATION THEORY 
A. Background 

There exist known important cases of critical phenomena, 
in which the infrared (IR) singularities one is concerned with 
originate from the divergence of a single length, such as the 
bulk correlation length . The most familiar example is 
static critical behavior that occurs at a usual bulk critical point 
when it is approached from the disordered phase. The crux of 
RG-improved perturbation theory is to exploit the RG to re- 
late the behavior of systems with large ^ to that of systems 
with £oo of order unity and then determine properties of the 
latter by means of appropriate perturbation theory techniques. 
Provided no other sources of IR singularities exist, the result- 
ing perturbation expansions will not be plagued by problems 
even though sophisticated resummation methods may well be 
required to obtain reliable results. 17 

In order to determine the behavior of properties at the 
critical point, information about the behavior of the static 
bulk analogs of the above scaling functions H± in the limit 
x /£oo - *■ is required. In the simple case considered so far, 
this can be inferred from the requirement that a temperature 
independent nonzero limit results for the property considered. 
Furthermore, the form of corrections to the leading asymp- 
totic behavior of the pair correlation function for |a;|/^oo <C 1 
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can be determined via the short-distance expansion (see, e.g., 
Refs. 68, 69, and their references). 



B. Breakdown of e expansion at T Cj00 and construction of 
effective action 



However, even in the case of static bulk critical behavior, 
additional sources of IR singularities may exist. For exam- 
ple, in systems exhibiting spontaneous breakdown of contin- 
uous symmetries, IR singularities associated with Goldstone 
modes appear on the coexistence curve. 70 Other examples 
are provided by crossover phenomena such as the crossover 
away from a bicritical point. Here the crossover from one 
type of critical behavior to another leads to singularities in the 
corresponding crossover scaling functions. 69 ' 71 To ensure that 
these IR singularities are properly taken into account in RG- 
improved perturbation theory is a nontrivial task. For some 
cases, specially designed RG schemes exists that allow one to 
correctly build in the asymptotic behavior at the stable fixed 
point to which the crossover occurs. 68 ' 72 Within the frame- 
work of dimensionality expansions, this usually works when 
the additional IR singularities are properly handled by sim- 
ple theories such as the random phase approximation 70,72 or 
are accessible to an e = d* — d expansion about the same 
upper critical dimension tf*. 69,71 Much more challenging are 
problems involving two distinct kinds of nontrivial critical be- 
havior with different upper critical dimensions. 

Unfortunately, the study of critical behavior in films of fi- 
nite thickness belongs to the latter class of hard problems. The 
reason is that a full treatment would involve a proper analy- 
sis of dimensional crossover. Suppose the film undergoes for 
finite L and given bulk dimension d a sharp phase transition 
at a temperature T c .l. If all interactions are ferromagnetic, 
this temperature must satisfy T c ^ < by Griffiths-Kelly- 
Sherman inequalities. 73 ' 74 On lowering the temperature from 
an initial temperature T > T c oo , one will therefore first ob- 
serve d-dimensional critical behavior as L,£oo — > oo with 
L/£,oo 3> 1, which will cross over to (d — l)-dimensional 
critical behavior as the length £|| on which correlations along 
the film decay becomes much larger than L. 14 Since d- and 
(d — 1) -dimensional critical behaviors involve different up- 
per critical dimensions (d* = 4 and d* = 5, respectively), 
one cannot handle both of them by the same dimensionality 
expansion. In those cases where no sharp phase transition is 
possible for finite L, a rounded transition will occur at T c l. 
Even then, the e expansion must not be expected to correctly 
capture the behavior for T ~ Tc.l- Because of these and 
technical difficulties, previous investigations of static finite- 
size effects in films based on the e expansion 49-5 1,75-78 have 
focused on the case T > T c oo , even though RG-improved 
mean-field results for T < T c exist for (d = 3)-dimensional 
systems with bulk critical and tricritical points. 79-84 

For simplicity, we shall restrict ourselves here also to the 
case T > T Ci00 . We begin with a discussion of the breakdown 
of the e expansion at T c oo for periodic and sp-sp boundary 
conditions. Since the latter is the only case of free bound- 
ary conditions explicitly considered henceforth, we shall no 
longer keep track of the enhancement variables, fixing c\ and 
C2 at zero whenever we deal with this boundary condition. 



We shall use dimensional regularization. In a perturbative 
approach, the critical values f c and c sp then vanish. 7,66 ' 67 In 
any case, the response and correlation propagators for sp-sp 
boundary conditions satisfy Neumann boundary conditions on 
both planes (p = N-N) at zero-loop order. 

The eigenvalues kf n of —d\ for periodic (p = per) and 
N-N boundary conditions are given by i;™ = 2-KmjL 

with m S Z and km ^ = irm/L with m = 0,1,..., oo, 
respectively. 49 In both cases, one has the eigenvalue fcg = 0. 
As can be seen from the mode decompositions (3.4) and (3.5), 
the associated ko = modes give IR singular contributions to 
the free propagators Rl and Cl at the bulk critical temper- 
ature (f = 0) when ui = p = 0. This is an artifact of the 
zero-loop approximation, which predicts that a sharp phase 
transition occurs for finite L precisely at T Ci00 , yielding no 
shift T c .l — T Ci00 . For bulk and semi-infinite systems the con- 
tributions from these zero modes are negligible. However, for 
films of finite thickness L, we must be more careful. 

It is known from Refs. 50 and 51 that beyond zero-loop 
order, the k = component of the order parameter becomes 



critical at the shifted values Sf = —5t 
given by 



(per) 



and —St 



(sp-sp) 



■o(pcr) 



2-c ^(sp-sp) 



n + 2 T(l -e/2)C(2-e) 
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27r 2- e /2 L 2-e 



(4.1) 

where ((s) denotes Riemann's zeta function. In order to ob- 
tain a well-defined RG-improved perturbation theory at T Cj00 , 
we generalize the approach of Refs. 50 and 51 to dynamics. 
Writing 

<j>(y,z,t) = L- l l 2 V { yi t)+^{y,z,t) 
4>(y t z,t) = J L- 1 / 2 v(y,t)+^(y,z,f) (4.2) 



with 



L »L 

iP(y,z,t)dz= / $(y,z,t)dz = 0, (4.3) 



we split the fields 4>(x,t) and <p(x,t) into their zero-mode 
components ip(y, t) and tp(y, t) and their k ^ orthogonal 
complements if>(x,t) and tp(x,t). Upon integrating out the 
latter fields, we define an effective action Jes\(p, tp] for the 
zero-mode components by 



J ef f[<p,ip] = -mTr^e 



-J[j>+L- 1/2 <p,il>+L- 1/2 <p] 



jM^^-ln/e-^'^). .(4.4) 



where 



(4.5) 
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Here 



-ld-l„ 



6L 



ip ]ip + cp — Xcp 



V 2 



(4.6) 



is the zero-loop contribution to J^f [ip, cp}. The interaction 
part J int is given by 



Xu . 
— / dt 
6L 



i) 2 (<p ■ <p) + 2(<p ■ i/>)(<p ■ if>) 

+ <p 2 (-0--0) + 2(-0- ¥ j)( V -V) 
+ L- 1 ' 2 ^ ■ if) yj 2 + iT 1/2 (y • V) ^ 

Computing Jeff in a loop expansion and writing 

Jeff [<p, <f\ = Jjff 1 [<P, <f] + Jiff' [<f, If] + ■ ■ ■ , 

one easily derives the one-loop contribution 



(4.7) 



(4.8) 



■ Tr In 



Xu 



G 



(p) 

L,ip 



B 



(4.9) 



(cf. Ref. 20), where Gjfl denotes the free propagator in the 
k 7^ subspace for the respective boundary conditions p = 
per and sp-sp. Let Pq = |fco)(^o| be the projector onto this 
k = space and Qq = 1 — Pq. Then we may write 



Gl,4> — QqGlQo 




RL,,ij(yi,ti; y 2l h) C L ^{y x , h; y 2 , t 2 ) 



® (<W) , 

(4.10) 



where Rl,^> and Cl4> are given by the analogs of Eqs. (3.4) 
and (3.5) one obtains by restricting their summations over m 
to values m/0. 

The operator B is block diagonal in yt space. Introducing 
the familiar symmetric tensor 



S a /3~f5 

we have 



1 



(4.11) 



B = 8( yi - y 2 ) 5(h - i 2 ) 







(4.12) 



<p=<p(yi,ti) 
<p=<p{yiM) 



The graphs (i)-(ii) depicted in Fig. 1 are the contributions 
implied by that involve two and four fields, respectively. 
Unlike graph (i), which is local in y and t, graphs (ii)-(iv) 
depend on the separations in y and t between their two ver- 
tices (marked as black dots) and thus are nonlocal. Similar 



y \ 
( 1 



(i) 



:* •'. 

s ✓ 

(ill) 



<p - 



(ii) 



(iv) 



FIG. 1: Graphs contributing to minus the effective action (4.8). Lines 
with and without an arrow represent the response and correlation 
propagators Rl and Cl, respectively. Red dotted lines indicate 
k — components; blue broken lines k 7^ components (color 
online). For further explanations, see main text. 



nonlocal contributions involving 6, 8, . . . fields are contained 
in jj^ [tp, tp]. At one-loop order no modification of the On- 
sager coefficient A appears. However, the two-loop graph (iv) 
produces a nonlocal modification of it. 
Let us use the notations 



(k;k) , - r 

75 1 ,...,a fe (yi>*l 



■ ■ • \Vk,tk) 

s' k+k j eS 



(4.13) 



and 



write 



^ai(y<ii,ii) • ■ -Sip ak (yk,tk) 

for the effective vertices of J e ff, 
7 s t (l/i j t±; . . . 1/fe, tk) for the analogously defined static 
effective k -point vertex functions (which were denoted by 
ry(k) j n R e f 51) gy analogy with the bulk case, the vertex 
functions 7^ 1;1 ^(. .. ;uj) and 7^ 1;3 ^(... ; {u>i}) must satisfy 
the relations 



and 



7 (1;1) (-.- ;c = 0) = A 7 i 2) (...) 
7^(...;{^ = 0}) = A 7 i t 4) (...), 



(4.14) 



(4.15) 



where the ellipses stand for the positions y t or momenta p^. 
The contribution of graph (i) to 7( 1;1 )(. . . ;u>)/X is indepen- 
dent of ui and hence agrees with that of its static analog 



1 1 



to the static two-point vertex function 7 S \ ' . It 
gives an L-dependent shift of the temperature variable, chang- 
ing it to 



where the integrals 



I^(L;r) 



dz \.f^\z)\ 2 

L P 2 + k^+T 



for the two boundary conditions in question, p 
sp-sp, are given by 



(4.16) 

(4.17) 
per and 



7 (per) {L . f } = f (d _ 3)/2 _ ^ f(d _ 2)/2 + 



2Q da {rL 2 ) 



(4.18) 



9 



and 



/f p - sp) (i;f) = / 1 (pcr) (2L;f) , 



(4.19) 



respectively. Here we started to use the notational conventions 
summarized in Appendix A for momentum integrals such as 
f . Further, 



A d 



2N, 



4-d 



± = -(4^)- d / 2 T(l-d/2) 



(4.20) 



while Qd,2{r) is a special one of the functions defined by 



85,86 



QdA r ) = 2 



E 



dk 

2^ 



(p2 + fc 2) CT /2-l 



P 



fc 2 



(4.21) 

Using the definition (4.16), the result implied by 
Eqs. (3.20)-(3.24) of Ref. 51 becomes 

7< 1;1 >(p,w)/A=p 2 -iw/A + 4 p) (f)+0(u 2 ) • (4.22) 



with such functions (namely, tpip -pvopeigcitoTs). Hence, the 
RG equations (4.24) carry over to the case of sp-sp boundary 
conditions. 

They can be exploited in a standard fashion to obtain the 
scaling forms of the vertex functions ^<( fc > fc ) ({p,, u^}) for both 
types of boundary conditions p = per and p = sp-sp. One 
obtains 

7 f ;fc >({p i)Wi };T,L)«A/A- 2 

x (^oo)-^-^XW ({p^; Wi/fi e }; L/^oo) (4.25) 



with 



- ^— -(fc + fc-2) 



7 " 2 
?? ; = (A + *)»j/2+(*-l)(a-2)/2, (4.26) 

and the characteristic bulk frequency 

fi c = A (m£oo)~ 3 • (4.27) 



For a discussion of the properties of the functions Qd.a the 
reader is referred to Appendix D of Ref. 51, where also plots 
of Qi,2{r) an d Q§.2{r) are displayed (see Fig. 9 of this ref- 
erence). It is known from previous works 49 ' 51 that the above 
results yield the shifts (4.1) for f = 0. This is easily verified 
utilizing the fact that 

lim Q da {r)/r = tt^" 5 ^ 2 T[(3 - d)/2] C(3 - d)/A 

r— >0 

= 7r~ d/2 r(d/2- l)C(d-2)/4 (4.23) 

(when d 7^ 1,3,5) according to equation (B35) of Ref. 87 and 
(C9)or(A15)ofRef. 51. 

It is equally easy to verify that the contribution of graph 
(ii) to 7C 1 ' 3 ) is in conformity with Eq. (4.15). Upon 
expressing in terms of 9tCx via the fluctuation- 

dissipation relation (2.11), one arrives at the desired result 
f™dtR L (t)C L (t) = (2A)- 1 [C L (t = 0)] 2 . 

Let us also note that the RG equations (3.15) for the func- 



tions W 



(N;N) 
R 



imply analogous ones for the renormalized ver- 
tex functions 7# . In the case of periodic boundary con- 
ditions, these vertex functions are evidently multiplicatively 

renormalizable. One has "f^'^ = Z 1 ^ 2 Z^ 2 "f^' k \ which 
leads to the RG equations 



2 



V4> 



(k;k) 

Tr 



0. 



(4.24) 



For p = sp-sp things are somewhat more subtle. It is 
known that r' 1 ' 3 ) has uv singularities corresponding to coun- 
terterms located on 03 that are proportional to <fi ■ d n <f> and 
<f> ■ d n cj). Therefore, it is not multiplicatively renormaliz- 
able and does not satisfy a homogeneous RG equation. 7 ' 10 ' 66 ' 67 
However, when integrated with smooth background fields sat- 
isfying the boundary conditions (2.12) with b\ — c 2 = 0, the 
boundary singularities of the mentioned form give no contri- 
bution. The construction of Jes[ip, ip] involves integrations 



V. SCALING FUNCTIONS OF FINITE-SIZE 
SUSCEPTIBILITIES 

We are now ready to turn to the computation of scaling 
functions by means of the small-e expansion. We first con- 
sider the scaling function of the dynamical response function 
\(tp ■ ip) at the bulk critical point. 



A. Dynamical finite-size response function 

Noting that the 0(n) symmetry is unbroken when T > 
T c .oo, we set a = (3 = 1 in the response and vertex func- 
tions and 7„0 to obtain 

-{V"P) = T t dz C dz'W$ 1] . (5.1) 
n L Jo Jo 

The pw-transform of this quantity gives us the finite-size sus- 
ceptibility 

Xl(p,u) = J Q ~J~ J dz ' 2 J dd l y^ J dtl2 



<5(^i(a;i,fi)) 
5hi(x 2 ,t 2 ) 



> j(oitl2-p-3/12) 



(5.2) 



The perturbation series of its inverse to order o(u 2 ) can be 
written as 



r (a), m 1 ~~ 9 o(e>) n + 2 u 



X e i) (0|f)+^). (5-3) 



where 



G^ 1 )(0|f) = -^_ 1 f( d - 3 )/ 2 



(5.4) 
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is the static bulk propagator in d — 1 dimensions 

G^- 1} (y|T) = / ( P 2 + f)-V™ 

J r> 



{r/y' 



Jd-3)/4 



'( 27r )(d-l)/2 A '(d-3)/2 (2/ 



/V^) , (5.5) 



taken at zero separation y. 

From the RG equations (3.15) one concludes that the renor- 
malized function Xlr nas as y m Ptotic scaling form 

X ^ R (p, u>] t) « M -"L 2 -" bCp) (pL, u/u) L , L/U) , (5.6) 

where we have introduced the characteristic finite-size fre- 
quency 



3 



(5.7) 



Using the above result (5.3), one can verify that this scaling 
form complies with the small-e expansion. At the order of our 
present calculation, the renormalization factors Z$, Z\, and 
Z u can all be substituted by unity, but we need 7 



Z T = 1 



3e 



u + 0(u 2 ) 



(5.8) 



to first order in u. The pole term cx fii/e of fx gets cancelled 
when f and u are expressed in terms of renormalized quanti- 
ties and the contribution from the counterterm oc Z T Z$ — 1 is 
added. The associated renormalized shifts 



n = n 

are uv finite; we obtain 

^(per) _ r (sp-sp) + 0(m2) 



(5.9) 



'L/2 



G 



2 u 



r e/2 



-4, 



d-1 



6/2 

fO(« 2 



{fi 2 L 2 T) d / 2 N d 
Combining the above results then gives 



(5.10) 



where 



-ILOfJb 

X 



^ + 0{u 2 ) , (5.11) 



» 



1 



+ 0{u 2 ) 



2 u^v 



3-e 



-(1+0/2 



6 fJ,LN 4 



In"'] 



} 

(5.12) 



are the inverse static finite-size susceptibilities. 

Note that by keeping the 0(u) terms of in the contri- 
butions ocu^]' 1 - 6 )/ 2 , we have included contributions of 

order u 2 in these terms. While these results for differ 
from those of Ref. 5 1 through precisely such terms, one eas- 
ily checks that they reduce to Eqs. (4.31) and (4.32) of this 
reference when expanded to first order in u. 



Our rationale for keeping the 0(u) contribution to the shift 
in the one-loop term associated with the zero mode should be 
clear. It prevents the propagator ^ that it involves from 
becoming massless at bulk criticality when L < oc. To com- 
pute the scaling functions E^ p \ we must compute x^f* a * t ne 
IR stable root u* = 3e/(n + 8) + 0(e 2 ) of the beta function 
j3 u , where we can exploit the facts that both 3 — 2 and r\ are of 
order e 2 . When r > 0, the shifted temperature variable f^ 
is of zeroth order in u. If we insert it into the contribution 
cx u[T^ ) ]( e - 1 )< /2 > the 0(e) term of r£ p) |„ =u . will therefore 
produce an 0(e 2 ) contribution. However, at criticality r = 0, 

f^ is linear in u and hence of order u* = 0(e) when eval- 
uated at u*. Accordingly, a contribution of order e 3/2 results 
from this source. 

Proceeding as indicated above, one arrives at the results 

^(p,ui,L)]- 1 = p 2 -iw + X^\L){l 

n + 2 



2vre^±| [^)(L)]- 1/2 } (5.13) 



with 

X(P er )(L) = 4l( sp - sp '(L/2) 

= L 2 + e!i ±2 16 7 r 2 Q 4 ,2(L 2 ) + 2 7 rL3 
n + 8 L z 

The latter functions have the asymptotic behaviors 

X^(l) « L 2 [1 + 0(L- X )] (5.15) 

L— *co 

and 



2tt 2 n + 2 



l^o 3a 2 n + 



(5.16) 



where a p is given by 




(5.17) 



Setting p = w = in Eq. (5.13) gives us the scaling func- 
tions of the static finite-size susceptibilities. One easily veri- 
fies that our results for L 2 /-^ (0,0, L) are in conformity with 
the 0(e) expressions (4.46) and (4.47) of Ref. 51. 

The limiting behavior (5.16) implies that the scaling func- 
tions (5.13) become 



— V6a p 
+ (e 3/2 ) 



3a2 n + 8 
+ 2\i/2 



.3/2 



(5.18) 



at L = 0, which explicitly shows the announced 0(e 3/2 ) con- 
tributions. Before we embark on a discussion of this result, let 
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us note that the scaling functions of the correlation function 
C^M = [ L ^- I" ' dz 2 [ d d ^y l2 I dt 12 



jo nL vo J J 

(5.19) 

need no separate calculations since the fluctuation-dissipation 
relation (2.1 1) in conjunction with Eq. (5.6) yields 



w 



(5.20) 



Had we computed the expansions (5.18) merely to linear 
order in e, one might be tempted to think that even the direct 
evaluation of these truncated series at d = 3 would give ac- 
ceptable estimates. The 0(e 3/2 ) terms reveal that there is no 
reason for such optimism: their signs differ from those of the 
contributions linear in e. Furthermore, their coefficients are 
quite large. Thus, evaluating the series expansion (5.18) to or- 
der e 3 /2 at e = 1 yields (unphysical) negative values even for 
n = 1, albeit the results remain positive when e is sufficiently 
small. This shows that more sophisticated extrapolation tech- 
niques are needed to obtain reliable estimates for these quan- 
tities at d = 3. 

An important check of the results for periodic boundary 
conditions can be made by comparing them with the exact 
large-n results obtained from the solution of the mean spher- 
ical model with short range interactions. 16 ' 86 ' 87 As we will 
show now, the large-n limit of the small-e expansion (5.18) 
for p = per can be recovered by solving the self-consistent 
equation that the finite-size susceptibility obeys in the limit 
n — » oo for d = 4 — e in an iterative fashion. 

There are two reasons why a similar comparison with large- 
n results cannot be made here in the case of sp-sp boundary 
conditions. First, translation invariance perpendicular to the 
boundary planes is broken for these boundary conditions, just 
as it generally is for free boundary conditions. The large- 
n limit therefore corresponds to a modified spherical lattice 
model involving separate constraints 88 ' 89 for the averages of 
Sigiayerz &f in each layer z, where Si is a spin variable on 
site i. Second, when d = 3, the thermal fluctuations should 
prevent the occurrence of film and surface phases with long- 
range order at temperatures T > in the continuous symme- 
try case n > 1 and hence for n — ► oo. Hence sp-sp boundary 
conditions cannot be realized when d = 3 and n = oo. 



B. Comparison with exact spherical-model results for periodic 
boundary conditions 

An important check of our results for periodic boundary 
conditions is their comparison with exact spherical-model re- 
sults. The static finite-size susceptibility of the spherical 
model on a <i-dimensional slab with periodic boundary con- 
ditions can be written as 



(SM) 

Xl 



where the scaling function X| M (L) is a solution [see, eg., 
equations (5.2) and (4.66) of Refs. 51 and 87, respectively] 



2Q dfl (Xf*)/Xf A = A ll [(Xf*) 



SM n (d-2)/2 d _ 2 i 



(5.22) 



For d = 3, the solution can be obtained in closed form; 87 ' 90 
it reads 



X| M (L) = 4arccscli 2 (2e 



-L/2\ 



4 In 2 



3 L/2 



(5.23) 



Its critical value, which corresponds to the amplitude 

l/S ( P er) (0,0,0),is 86 



X 



3 SM (0) = 41nV 



0.926259, g = 1+V5 )/2 , (5.24) 



where g is the golden mean. 

To determine the amplitude X^ M (0) in d = 4 — e di- 
mensions, we substitute the series expansion of the function 
Qd, 2 { r ) given in equation (C.9) ofRef. 51 to obtain the repre- 
sentation 



l 2 -A^ t r^' 2 + Ur) (5.25) 



with 
where 



fe(r) 



J2 bi j) rj 

3=0 



i-iy t[j + ( £ - 1)/2] 

j\ 2 2 J +1 7T 2 i+( € + 1 )/ 2 



3 



C(2j 



(5.26) 



(5.27) 



The coefficients ^44_ e and bi^ have simple poles at e = 
whose residues differ by the factor (—2), so that the pole term 

of Ai_ c r 1_£/2 cancels that of be r. The remaining coeffi- 
cients, ^3_ e and be with j =/= 1, are regular at e = 0. On the 
other hand, /l3_ c and be have simple poles at e = 1 whose 
residues are equal. 

These results imply that X^_ t (L) — L 2 must vanish linearly 
as e — > 0, so that the right-hand side of Eq. (5.22) approaches 
a nonzero limit. The 0(e) term of X4_ e (L) — L 2 was com- 
puted in Ref. 51 [see its Eq. (5.3)]. To determine the 0(e 3 ' 2 ) 
contribution to X4_ e (0), we set L = and solve Eq. (5.22) 
with the ansatz X4_ e (0) = C\t 
gives 

2tt 2 



C 3/2 e 3/2 + o(e 3/2 ). This 



L'/Xr(L/U) 



(5.21) 



Xl M M = =He - V6e 3 / 2 ) + o(e 3 / 2 ) . (5.28) 

Obviously, the result agrees with the expansion (5.18) of 
[HO^O^O)]" 1 forn = oo. 



C. Dynamical finite-size surface response function 

In Sec. V B we were concerned with the dynamic finite- 
size susceptibilities (5.2). These are integral quantities. As 
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an alternative, we here wish to consider the local susceptibil- with 



= \(<p^(y,t) ■ 4>^(y,t))/n and X [V P) 



ities X <r P) 

A (</> <Bl (y, t) ■ <p s ' 2 (y,t))/n describing the linear responses 
of the order-parameter density at one boundary plane to a 
magnetic field acting on the order-parameter density in the 
same and the complementary boundary planes, respectively. 
In the case of periodic boundary conditions, both quantities 
are the same — and by translational invariance along the 
z-direction — identical with the layer susceptibility Xzz = 
\(<p(y,z,t) ■ <p(y,z,t))/n for any layer z. We start with 
an investigation of the latter. For the sake of simplicity, we 
do not consider deviations from the bulk critical temperature 
here, setting r = 0. 



1. Layer susceptibility for periodic boundary conditions 

Our analysis of the finite-size susceptibility xl suggests 
that we should work with a dressed (^-propagator that ac- 
counts for the 0(u) shift f| per ' ) given in Eq. (4.1). Noting that 
the ko = 0-mode contribution to Cl{p — 0; z = 0; t = 0) 
vanishes when f = 0, one sees that a single two-loop graph 
suffices to capture all contributions to order e 3 / 2 in the renor- 
malized quantity Xzz.r- One obtains 



Xzz(p,u;L)/\ 



coth(k u L/2) 
2\n u] 



O < s m' < O 

z z 



= \]p 2 - ii 



(5.32) 



The susceptibility Xzz.b. must reduce to an L-independent 
bulk quantity as L — > oo. Therefore, the scaling function 
should have the limiting behavior 



"(per) 



(p,w) « p''" 1 *(wp~ 3 ) . (5.33) 

p — >oo 



Equation (5.33) shows that this is indeed the case, giving 

1 



*5?(wp 3 ) 



2yl — iwp~ 



+ 0(e 2 ), (5.34) 



which is the correct bulk result. 

As long as L is finite, there is no reason for Xzz.n to di- 
verge at the bulk critical temperature. Hence, we expect 
its inverse 1/xzz,r(p, 0; L, 0) to approach a nonzero limit. 
From Eq. (5.31) we can read off the small-e expansion of 

1 /E« r ^ (p, w) and compute its limit 



lim 



1 



p-0 ~i per) (p,0) 



2?r 2 n + 2 
~3~ n + 8 

+ (e 3 / 2 ) . 



e-A/6 



,3/2 



(5.35) 



0+0(u 2 ) , (5.29) 



O « • « 



where k u now denotes the quantity (3.10) with f set to zero. 
Here full black lines with and without arrows represent the 
free response and correlation propagators and Cl with 



-i (per) 



with 



f = 0, respectively. The dashed blue line indicates , 
f = 0. The dotted red line is the k = part of C^ pel ^ with 

f = <5fr p . Note that we do not include the shift 5rj pel ^ in 
the k = contributions to the external legs of the two graphs. 
Their inclusion would produce terms of order u 2 and hence 
lead to corrections of order e 2 . 

In Appendix C the graphs in Eq. (5.29) are computed. From 
the results given in Eqs. (C2) and (C3), the RG-improved per- 
turbation expansion of the renormalized dynamic layer sus- 
ceptibility Xzz.b. = Xzz/Z^ follows in a straightforward fash- 
ion. Evaluating it at the fixed-point value, one sees that it is in 
conformity with the scaling form predicted by the RG equa- 
tions (3.15), namely 



S(p^;£,t = o) 



^""S^CpL.w/wl). (5.30) 



where the scaling function has the small-e expansion 

coth(x/2) 7r 2 n + 2 x + smhx 



2x 



e-A/6 



12 n + 8 ^sinh 2 ^) 

o(e 3 / 2 ) (5.31) 



J3/2 



The result agrees at this order with the expansion (5.18) 
of the inverse static susceptibility [S' per ^(0, 0, 0)] _1 . Thus 
the inverse static zero-momentum layer susceptibility 
1/Xzz,fl(0, 0; L, 0) is indeed nonzero. Of course, the previ- 
ously discussed difficulties to obtain reliable d = 3 estimates 
from this expansion to the order e 3 / 2 apply here as well. 



2. Surface susceptibilities for sp-sp boundary conditions 



and 



The graphs of the surface susceptibilities Xn 

X^l sp ' corresponding to the ones displayed in Eq. (5.29) 
are computed in Appendix C 2. The results are gathered in 
Eqs. (C25)-(C27). Using them one can verify in a straightfor- 
ward manner that the renormalized susceptibilities Xn It = 

x[T Sp) Wi and - X^/Z^Zi ^ uv rmite'to the 

appropriate order in u. One obtains 



..(sp-sp) 
All B 



coth(K^L) n + 2 f coth(K w L) 



K u 6 ^ 

2InO*L)] + S(k u L) + R(k u L) 



l-2 lL 



y 



, 4 / n+2 



1/2-, 



0(e) 



(5.36) 
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and 



(sp-sp) 
XlL.R 



csch(K w L) n + 2 (csch(K w L) 
— u< [1 - 2^ E 

- 2 lnGuL)] + 2S{ 0) 2 (k u L) + R 2 {k u L) 



IT 

y 

0(e) 



, A f n + 2 
1 — - I - „ ,„ u 



L \2fi 2 L 2 



1/2 



(5.37) 



where k u 

?(0) 



yV 2 — iuin 2 /\, while 5(>r), R{x), ii,i(x), 

S^l{x), R2(x), and l\ t i(x) denote the functions specified 
by'Eqs. (C22), (C23), (C8), (C20), (C24), and (C9), respec- 
tively. The functions R(x) and R-^ix) are single integrals, 
which can be computed by numerical integration. 

According to the RG equations (3.15) these susceptibilities 
should have the asymptotic scaling behavior 



X ^J } (p, lj; L, r = 0) « /T^MF Xn ( P L, u/uo L ) , 

1L.R 1L 



(5.38) 

where ?/j ,p is a standard surface correlation exponent asso- 
ciated with the special transition whose e expansion rf^ = 

-e (n + 2)/(n + 8) + 0{e 2 ) is known 7 ' 66 ' 67 to order e 2 . Eval- 
uating the above results (5.36) and (5.37) at u* shows their 
consistency with these scaling forms and yields the small-e 
expansions 



X u (p,w) 



coth k n + 2 e 



x n + 8 2 
+ S(x) + R{x) 



o (1 - 2 7 £) 



coth; 



7T 

y 



1 - 2 6e 



2\!/2' 



and 



Xil{p,w) 



+ o(e 3 / 2 ) 



cschx n + 2 t 



(5.39) 



x n + 8 2 

7 (0) 



(1 - 2 7 i0 



csch>r 



+ 2^(x)+i? 2 (x) 

n + 2\i/2 



^ 2 

+ y 



1 - 2 6e 



V n + 8/ 



+ o(e 3 / 2 ) 



(5.40) 



where x is again given by Eq. (5.32). 

In the limit p — > oo, the function Xn must behave as 



X 11 (p,w) w p"* l *n(«;p-») 



(5.41) 



where ^ n must agree with the scaling function of XiT h °f 

the semi-infinite (L = oo) theory. This is indeed the case. The tibilities at T CiOC to order e 3/2 . We were also able to cast the 



limit liirip^oo p (p, w) can be computed in a straight- 



forward fashion, giving 



#11 (u) 



1 



1 



e n + 2 



[l-ln(4-4w)]}+0(e 2 ) 



v/1 - w I 2 n + i 

(5.42) 

An independent calculation for the semi-infinite theory yields 
precisely the same result. 

oo, the analog of the seal- 



Since Xil S r vanishes as L 



ing function ipu for xtl r must vanish. Checking the limit 
#ii(p, w) explicitly, we do in fact find that it 



lim p ^ 
is zero. 



1— Tf} 
> ''II 



Just as the layer susceptibility Xzz,R> the surface suscepti- 
bilities xf 1 p "Jf' ) and Xil'^r should not become critical at the 
bulk critical point and hence have finite p — > limits at 
t = co. To check this, we set u) = in Eqs. (5.39) and 
(5.40), determine the series expansions of [Xn(p,0)] _1 and 
[X 1L (p, 0)}- 1 to 0(e 3 / 2 ), and then take the limits p -> 0. The 
results are 



[Xii,«(o, o)]- 1 = [X 1LtR (0, 0)]"^^ 3 / 2 ) 



ir 2 n + 2 
y n + 8 



1 - 2 6e 



2\ 1 /2' 



n + 8 



o(e 3/2 ). (5.43) 



They indicate that an i-dependent shift of the multicritical 
special point occurs also along the c-direction, not only along 
the temperature direction. Note, however, that a special point 
at which a bulk disordered phase with long-range surface or- 
der meets with a bulk-disordered, surface-ordered phase and 
a bulk-ordered, surface-ordered phase exists for d = 3 di- 
mensional, semi-infinite spin systems only when n = 1. The 
analog of this multicritical point for finite L is the one where 
the transition temperature of the surface transition coincides 
with the transition temperature T Cj l of the film. 



VI. CONCLUSIONS 

In this article, we have studied model A in film geometry. 
We focused on the cases of boundary conditions for which the 
classical (Landau-van Hove) theory involves zero modes: pe- 
riodic boundary conditions and free boundary conditions with 
critically enhanced surface interactions. 

Major motivations were to determine the consequences of 
the zero mode and to reformulate RG-improved perturbation 
theory such that a well-defined small-e expansion results at the 
bulk critical temperature. Our reformulation of RG-improved 
perturbation theory builds on the strategy pursued in Refs. 50 
and Ref. 5 1 to investigate static properties such as the finite- 
size free energy and the Casimir force. In these papers it 
was shown that the conventional e-expansion breaks down and 
contributions involving half-integer powers e fc/2 with k > 3 
modulo powers of In e appear in the expansion of the critical 
Casimir amplitudes and other quantities. Here, we confirmed 
this breakdown and explicitly determined the small-e expan- 
sions of various dynamic finite-size, layer, and surface suscep- 
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e-expansion of the temperature-dependent scaling functions of 

the finite-size susceptibilities Xlr f° r both boundary condi- 
tions p = per and p = sp-sp in a form that reproduce their 
expansions to order e 3/2 as T — > T coo . 

Despite these successes, the theory suffers from severe lim- 
itations and is still in an unsatisfactory state. First of all, our 
results (5.35) and (5.43) for the inverse layer and surface sus- 
ceptibilities at T c oo suggest that reliable estimates for d = 3 
can hardly been obtained on the basis of such small-e expan- 
sion to order e 3/2 alone. Just as in the case of the Casimir 
amplitudes and the scaling function of the residual free en- 
ergy at small values of L/£oo, 51 the deviations of the simplest 
extrapolations obtained by setting e = 1 in the series truncated 
at the respective orders e and e 3 72 appear to oscillate in sign. 
It is conceivable that the situation is particularly bad in those 
cases where the classical theory involves zero modes. These 
zero modes provide a separate source of IR singularities, a 
problem one also encounters in thermal field theory. 91 It has 
been suggested to handle IR problems of this kind by a resum- 
mation of "foam" diagrams. In the case of the n-component 
bulk <j> 4 theory at finite temperature, 92 this procedure yields 
results for the pressure in conformity with the exact large-n 
result. Whether and to what extent RG-improved perturbation 
theory might be improved by combining it with resummations 
of this kind remains to be seen. 

Finally, let us emphasize that there is little reason to believe 
that the theory is in a much better state in those cases where 
the classical theory does not involve zero modes, notwith- 
standing the additional difficulties such modes cause. For free 
boundary conditions, one will generically encounter a zero 
mode in the classical theory at an L-dependent temperature 
T^l jtz T Cj00 . This indicates that in Landau theory the film 
becomes critical at this temperature and undergoes a transi- 
tion to an ordered low-temperature phase. The ordered phase 
and hence the transition to it may not survive the inclusion of 
thermal fluctuations when L < oo. This happens in the con- 
tinuous symmetry case n > 1 when d < 3, where it should 
be recalled that the d = 3 case with n = 2 is special in that 
the film has a low-temperature phase with quasi-long-range 
order. However, even in those cases where the film does have 
a transition to an ordered low-temperature phase for finite L 
(as it does when d = 3 and n = 1), one encounters two im- 
portant challenges that are beyond the scope of the presently 
available analytical RG approaches but any satisfactory full 
theory of dimensional crossover must be able to cope with: (i) 
to determine the location of the singularity of the residual free 
energy's scaling function corresponding to the transition tem- 
perature T Cj £ with acceptable accuracy, and (ii) to yield the 
correct IR singularities at this transition in conformity with 
the expected (d— l)-dimensional critical behavior of the film. 
The difficulty is that even the shift cannot normally be com- 
puted by perturbation theory but requires RG techniques to 
deal with the IR singularities. 93 The RG scheme employed 
here and in the work of Krech and Dietrich 49 may be appro- 
priate to go on scales of the order of ^ as long as < L. 
However, it is insufficient to integrate out degrees of freedom 
between L and the film's correlation length in an adequate 
fashion when £j, > L, and to correctly yield the IR singular- 



ities at the critical temperature T c ^ even when the boundary 
conditions do not involve zero modes at T Ci00 . 

We close with some comments on the universality of our 
results and finite-size scaling results in general. Chen and 
Dohm 56 ' 57 recently launched a discussion of the universality 
of finite-size scaling results and the validity of two-scale (and 
multiscale) factor universality. Let us consider their concerns 
in some detail. A first issue raised in Ref. 56 is that the use of 
a sharp large-momentum cutoff modifies the L dependence of 
the singular part of the finite-size free energy density of sys- 
tems of linear size L in a qualitative manner. This effect is 
unphysical and entirely due to the use of a sharp cutoff; that 
a sharp cutoff can produce unphysical effects has been known 
since the early days of Wilson's RG. 2 For systems of the kind 
considered by them — systems that are finite in all directions 
— the issue was discussed and clarified in Refs. 94 and 96. It 
needs no further discussion. 

A second point made in Refs. 56 is that long-range inter- 
actions which are irrelevant in the RG sense produce alge- 
braically decaying contributions to the singular part of the 
finite-size free energy density. Such interactions were previ- 
ously considered by Dantchev and co-workers, 94 - 95 who intro- 
duced the term "subleading long-range interaction" for them. 
Away from criticality, these contributions compete with the 
exponentially decaying ones one has for systems with purely 
short-range interactions and become dominant in the appro- 
priate region of temperature and large L. While such contri- 
butions (which are expected to be small on an absolute scale) 
still have to be clearly verified by experiments, they certainly 
are real. 

Chen and Dohm 56 interpreted their presence as signaling 
the breakdown of finite-size scaling. However, what is broken 
is just the simple version of finite-size scaling that involves a 
single length besides L, namely, the correlation length. Sub- 
leading long-range interactions give rise to at least one further 
length — the one associated with the corresponding irrelevant 
scaling field. 97 This must not be set to zero in order to retain 
the long-range tail in the regime where it dominates the expo- 
nentially decaying short-range contribution to the finite-size 
free-energy density. It may well be set to zero in the crit- 
ical regime £/£oo <C 1 because the subleading long-range 
interaction contributes there only a correction to the leading 
L dependence. Thus, subleading long-range interactions are 
intermediate between dangerous irrelevant and conventional 
irrelevant perturbations: they share with the former the prop- 
erty that they must not generally be set to zero. Unlike those 
(which would affect the leading critical behavior of some 
quantities), but similar to the latter, they give only corrections 
to the leading critical behavior. 

Note that the mechanism just described for subleading 
long-range interactions is neither specific to finite-size critical 
behavior nor new. A familiar analog known from the study 
of critical adsorption of fluids was discussed more than 25 
years ago by de Gennes. 98 Substrates ("walls") typically ex- 
ert one-body forces on the fluid that have besides short-range 
components algebraically decaying van-der-Waals tails. In a 
semi-infinite geometry bounded by a wall at z — and re- 
stricted to z > 0, the latter contribute effective wall-fluid 
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interactions of the form J d d x y J °° dz h w f(z) <j)(y, z) to the 
Hamiltonian, where h w f(z) behaves as A w fZ~ v as z — > oo 
(cf. Ref. 99 and Sec. 3.11 of Ref. 7). The long-range part is 
irrelevant in the RG sense provided the exponent v is larger 
than the magnetic RG eigenexponent yh = (d + 2 — rj)/v, 
which is the case for nonretarded and retarded van-der-Waals 
interactions in d = 3 dimensions (for which v = 3 and v = 4, 
respectively). It produces a long-range tail ~ z~ v to the de- 
viation 8m(z) = m(z) — m(oo) of the order-parameter den- 
sity m(z) = (4>(y,z)) from its bulk value m(oo). By set- 
ting the length ~ g^f~ Vh associated with the irrelevant scaling 
field g w f ~ A w f to zero, one would loose this algebraically 
decaying contribution. On the other hand, the leading temper- 
ature singularity ~ | r |-<>-' 3 ) of the excess order parameter 
J °° dz Sm(z) would remain the same since <? w f yields merely 
corrections to scaling for this quantity. The analogy with how 
subleading long-range interactions affect finite-size properties 
is obvious. Of course, the amplitude of the irrelevant scaling 
field <7„f is nonuniversal. Following the logic of Refs. 56 and 
57, one would have to call this a violation of scaling in semi- 
infinite systems, though it again just means that single-length 
scaling reaches its limits, failing to capture the asymptotic be- 
havior of certain quantities. 

To what extent would the inclusion of subleading long- 
range interactions alter the results of our analysis of dynamic 
finite-size critical behavior given above? While a detailed, 
quantitative analysis of their effects is beyond the scope of 
this paper, clear predictions can be made on general grounds 
and the basis of what is known from statics. Power laws de- 
scribing asymptotic dependences in L or r at T — T c oo will 
be modified by corrections to scaling involving the associated 
irrelevant scaling field. 97 For a subleading long-range pair in- 
teraction decaying ~ x -{d+a) ( w jflj a > 2 — if), the as- 
sociated correction-to-scaling exponent is ui a = a — 2 + r\ 
(see, e.g., Ref. 87 and its references). Thus these correc- 
tions should be down by factors L^^ a (or S,^") in compar- 
ison to the respective leading power laws. Correspondingly, 
scaling forms such as the one for the finite-size susceptibil- 
ity Xif w > r ) gi ven m Eq. (5.6) should obtain a correction 

oc L 2 ~ ri ~ UJ " Eff^{pL, (j/uji^L/^oo) at linear order in the ir- 
relevant long-range scaling field. 97 Furthermore, there exist 
quantities whose behaviors get qualitatively modified by sub- 
leading long-range interactions. This is typically the case (in 
certain regimes of L and r) for quantities that decay expo- 
nentially in the absence of long-range interactions. Obvious 
examples are zero-frequency response functions in position 
space at temperatures r > 0; these decay algebraically in the 
limit of large distances Xij = \xi — Xj\ — ► oo between two 
points when the pair interactions have a subleading long-range 
tail. 

Of course, a proper investigation of the effects of sublead- 
ing long-range interactions should also allow for irrelevant 
surface-related scaling fields (such as pair interactions local- 
ized on the boundary that decay algebraically as a function of 
the separation yt — yj along the boundary planes and pair in- 
teractions in the interior of the sample that decay as a power 
of the distance from the boundary planes). This is beyond the 



scope of our present work. 

We conclude by turning to a third source of universality vi- 
olations, discussed extensively in Ref. 57: the effects of weak 
anisotropy. To keep things as simple as possible, it will be 
convenient to discuss the issue first in the context of static 
bulk critical behavior. The characteristic property of systems 
exhibiting weakly anisotropic bulk critical behavior is that the 
correlation lengths describing the decay of correlations along 
arbitrary directions diverge ~ | | — ^ with one and the same 
critical exponent v, but the shape of the correlation region is 
ellipsoidal rather than spherical. This means that the square 
gradient term of the Hamiltonian (2.1) in general takes the 
form 

\ f d d xB kl (d k cf>) ■ drf (6.1) 

in Cartesian coordinates, where dk4> = d<p/dx k are par- 
tial derivatives with respect to these Cartesian coordinates 
x k , k = 1 , . . . , d, and Einstein's summation convention is 
used. The matrix B = (B kl ) is symmetric and positive defi- 
nite. Hence its inverse exists and defines a metric tensor B k i ■ 
Accordingly the modified square gradient term (6.1) can be 
viewed as the scalar product of the gradient operator with it- 
self in this metric. 

An evident consequence of the choice (6. 1) of the modified 
square gradient term is that a corresponding replacement 

V • V -> l) k B kl l) l (6.2) 

must be made for the second-order derivative operator in the 
dynamic action (2.4). 

In cases where the underlying microscopic system whose 
critical behavior one is concerned with has cubic or or- 
thorhombic lattice symmetry, this matrix B is proportional to 
the unity matrix 1 or at least diagonal, but for monoclinic and 
triclinic lattices it is generally nondiagonal. It can be trans- 
formed to 1 by combining an orthogonal transformation O 
with O 1 ■ B ■ O = diag(6 1 , . . . , b d ) = b to principal axes 
with a rescaling of coordinates. Let us make the coordinate 
transformation 

if : x = {x k ) ^ x' = (x' k = tp k (x\...,x d )) 

= b- 1 ' 2 O 1 x (6.3) 

and introduce the transformed quantities 





= B^ 4 4>(x,t) , 


4/(x',t) 


= £? 1/4 <fr(x, t) , 


t' 


= T , 


u' 




h'(x',t) 


= B^ A h{x,t) , 


A' 


= A, 


3? 





where B = det B. 

Consider the analogs of the Hamiltonian (2.1) and dy- 
namic action (2.4) with the modified square gradient terms 
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Eqs. (6.1) and (6.2), respectively. Using the definitions (6.4), 
we can express their contributions involving the volume in- 
tegrals d d x in terms of transformed ("primed") quantities. 
One easily verifies that the resulting expressions are identical 
with the original ones given in Eqs. (2.1) and (2.4), up to the 
replacement of unprimed by primed quantities. In particular, 
the primed square gradient term takes the standard form 



2 ./a 



^ dx' k dx' k 

k=l 



a r 



(6.5) 



All those transformations given in Eq. (6.4) that refer to 
static bulk quantities and the mapping of the finite-size re- 
gion 23 are consistent with those of Ref. 57 (whose matrix A 
corresponds to our B). The remaining ones are required for 
dynamics. To cope with free boundaries, we should also deter- 
mine which boundary contributions to the transformed Hamil- 
tonian and dynamic action result from the boundary integrals 
I<B d d ~ 1 y ■ ■ ■ of Eqs. (2.1) and (2.4), respectively. Before we 
do this, let us briefly discuss what it means for bulk critical be- 
havior that the Hamiltonian describing the asymptotic critical 
behavior of weakly anisotropic systems can be mapped onto a 
primed one with coefficient matrix B' = 1. 

There is no question that this mapping to primed variables 
involves nonuniversal parameters. It is also clear that for non- 
diagonal B kl , part of the nonuniversality resides in the direc- 
tions of the principal axes, as emphasized by Dohm. 57 Does 
this mean that fundamental concepts of the modern theory of 
critical phenomena such as the notion of universality classes 
for static bulk critical behavior must be questioned, crucially 
modified or even given up? We see no reasons for such a 
conclusion. In our view, the very existence of the above map- 
ping to simple minimal models such as the conventional 4 
theory is a clear signature of universality because it ensures 
that the critical properties of weakly anisotropic system can 
be expressed in terms of the universal properties of the latter. 
(Inasmuch as the explicit results of Ref. 57 are concerned, we 
are not aware of discrepancies; yet, our view of the situation 
may not fully be shared by its author.) 

The nonuniversal geometric dependences contained in the 
metric should come as no surprise. An essential element of 
the modern theory of critical phenomena is a mapping of mi- 
croscopic models onto conceptually simple continuum mod- 
els such as the (f> A theory whose critical fixed points describe 
the respective universality classes of static bulk critical behav- 
ior. Such mappings always involve nonuniversal parameters. 
Two familiar examples of such parameters are the location of 
the critical point and the slope of the coexistence curve; their 
nonuniversality shows up in the dependence of the two rele- 
vant scaling fields g T and on T— T Cj00 and the deviation of 
the magnetic field or chemical potential /i from their critical 
values H c oo = or ft CiCO ,' where it should be recalled that in 
the case of fluids, both the thermal and magnetic scaling fields 
are nontrivial linear combinations of T — T c oo and fi — fi C)00 , 
up to nonlinear contributions. 

The nonuniversal geometric dependences contained in the 
metric are of an analogous kind. Absorbing them through the 
choice of properly defined transformed quantities is similar to 



the adsorption of other nonuniversal properties such as the lo- 
cation of the critical point and the slope of coexistence curve 
via appropriately chosen scaling fields. Of course, in any com- 
parison of predictions of the theory with experimental results 
or Monte Carlo simulations for systems with weak anisotropy, 
the nonuniversal geometry associated with the metric must be 
taken into account since it enters the way lattice observables 
depend on the order parameter. 

In most studies of critical behavior either standard square 
gradient terms with B = 1 are chosen from the outset or else 
it is tacitly assumed that the above transformation to primed 
variables has been made. This is done, in particular, when 
two-scale-factor is defined and discussed. In Ref. 57 it is 
emphasized that two-scale-factor universality is broken un- 
less B = 1. Formally, this is correct since nonuniversal pa- 
rameters that cannot be absorbed in the nonuniversal ampli- 
tudes of the two relevant scaling fields g T and g^ are involved. 
However, we believe it is natural and more reasonable to de- 
fine two-scale-factor universality only after the transforma- 
tion to primed variables has been made. How the two-scale- 
factor universality of the corresponding conventional <f)' 4 the- 
ory manifests itself in the original one with non-Euclidean 
metric follows from the relation between these two theories 
(as is worked out in detail for the case of static critical behav- 
ior in Ref. 57). 

These considerations generalize in a straightforward fash- 
ion to the case of dynamic bulk critical behavior described by 
model A with B ^ 1. To analyze the corresponding asymp- 
totic dynamic critical behavior, a dynamic scaling field (asso- 
ciated with the Onsager coefficient A) must be considered in 
addition to g T and gn- Hence a third nonuniversal scale factor 
must be fixed in the corresponding transformed theory with 
Euclidean metric. 

A first obvious, though important, new feature one encoun- 
ters when extending these considerations to finite-size systems 
is that the integration region 23 and its boundary 923 transform 
under the map <p. (Note that momenta would transform with 
the inverse map v?" 1 , so integration regions in momentum 
space and hence momentum cutoffs would transform as well.) 
For general matrices B this is a shear transformation; cubes 
of finite linear dimension get transformed into parallelepipeds. 
Therefore, finite-size systems of a given (say, cubical) shape 
that involve the generalized square gradient terms (6.1) and 
(6.2) should not be compared to their primed analogs of the 
same but of a different (non-cubical) shape. 

According to the phenomenological theory of finite-size 
scaling, the finite-size critical behavior of a given microscopic 
system does not only depend on those properties that deter- 
mine the corresponding bulk universality class but also on 
gross finite-size properties such as shape and boundary con- 
ditions. Hence, each universality class for static bulk critical 
behavior generally splits up into several universality classes 
for finite-size critical behavior. This is analogous to the split- 
ting up of static bulk universality classes into several distinct 
universality classes for dynamic bulk critical behavior and 
into those for static boundary critical behavior. 7 - 8 The upshot 
is that two finite-size systems with the same volume region 23 
(and hence shape) whose large-length scale descriptions in- 
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volve generalized and standard square gradient terms, respec- 
tively, may represent distinct finite-size universality classes 
even when the same kinds of boundary conditions (e.g., pe- 
riodic boundary conditions) are chosen for both of them on 
the level of lattice models. 

Square gradient terms with nondiagonal B give rise to im- 
portant modifications of the finite-size critical behavior of sys- 
tems that are finite in all directions. This is discussed in detail 
for the case of static critical phenomena in Ref. 57. In the 
case of our slab geometry, the image 2J' of the slab of infi- 
nite lateral extension and thickness L under ip is again a slab 
whose thickness L' generally differs from L. However, we 
must also clarify how the boundary conditions are affected by 
generalized gradient terms and what happens to them under 
the mapping to the primed system. Considering a Hamilto- 
nian and dynamic action that agree with those specified by 
Eqs. (2.1) and (2.4) except for the replacements of (V</>) 2 /2 

and V • V by the generalized gradient term (6.1) and the oper- 
ator (6.2), respectively, one finds that the boundary conditions 
now become 

n T ■ B ■ V0(x, t) = cj <fi(x, t) , x G OS, , 

n T ■ B ■ \7<f>(x, t) = cj (j>(x, t), x e QSj , (6.6) 

where n T is the row vector transposed to the column vector 
n. 

Upon introducing the vector / = // = B ■ n, one recog- 
nizes the derivatives on the left-hand sides of these equations 
as directional derivatives df<p = fd^<p. Recall that in our 
case n = ±e z on $81,2- For general B, the vector / is not 
parallel to n. The meaning of these boundary conditions can 
be understood as follows. Suppose we extrapolate the fields 
4> and 4> from a point x<&. on boundary plane 58 j along the 
direction — / in a linear fashion. Then these extrapolations 
vanish when the coordinate differences (x — xr Sj ) ■ f take the 
values ±//cj for j = 1 and 2, respectively. Thus, the modi- 
fied square gradient term in general changes both the direction 
of the extrapolation and the associated extrapolation length. 
Further, the vector / is normal to the surface 23 in the met- 
ric (Bki). To see this, let us represent vectors v and w in the 
canonical basis {ej} ofM. d as / = pe.j and denote the scalar 
product in the metric (Bki) as B(v, w) = Bkiv h w l . The dual 
basis {/^ fe '}, which satisfies B{f^ h \ei) = 5 k i, is given by 
f( k ) = B kl ei. For points on Q3i.2, the vector / is nothing but 
±/W, and hence orthogonal to the tangent vectors e^d in 
the metric (Bki). In fact, the directional derivative df on the 
left-hand sides of Eq. (6.6) corresponds to a normal derivative 
in this metric. Transforming to primed variables, gives the 
Robin boundary conditions d' n 4>' = c'j4>' for the fields <$' and 

4>' with c'j = Cj/ /, in conformity with our above results. 

Since the enhancements variables Cj and c'j of the systems 
with generalized and standard square gradients are propor- 
tional to each other, the fixed points of the transformed sys- 
tem with c'j = 0, ±00 map onto the respective fixed points 
of the unprimed system; the effects of the generalized square 
gradient terms implied by weak anisotropy normally may be 
expected to be of a purely geometrical kind. They should be 
particularly important and interesting when B is nondiagonal. 



One class of systems deserving detailed studies consists 
of binary alloys. This is because their description gener- 
ally involves, besides the order parameter, further, so-called 
secondary densities (nonordering densities). Studies of the 
static boundary critical behavior of body centered cubic bi- 
nary alloys have revealed that careful investigations of the 
coupling of the order parameter to these secondary densities 
and the symmetry reduction caused by the presence of bound- 
ary planes may be necessary to determine which surface uni- 
versality class applies. In fact, depending on the orientation of 
the surface plane relative to the crystal axes, distinct surface 
universality classes may be realized. 101 ' 102 Analogous studies 
have yet to be performed for binary alloys with less symmetric 
(e.g., monoclinic and triclinic) crystal structures, which would 
yield nondiagonal metrical coefficients B kl . 

We close with a brief discussion of an elementary exam- 
ple of a slab exhibiting weakly anisotropic critical behavior. 
Consider a nearest-neighbor (NN) lattice O(n) spin model 
that is restricted to the layers z = 0, 1, . . . , L of the sim- 
ple cubic lattice Z d . To introduce weak anisotropy, we as- 
sume that all NN bonds perpendicular to the bottom and top 
layers z = and z = L have strength J±, while those 
along the top, bottom, and remaining layers have different 
strengths Jim, J\\ 2> J\\> respectively. Mapping this micro- 
scopic model onto a continuum model gives squared gradient 
terms with B = diag(6 y , . . . , b y , b z ), where b z /b v = Jj_/Jy. 
Hence the ratio of the corresponding bulk correlation lengths 
£oo,\\ an d Coo, _i_ (defined via second moments of the respec- 
tive displacements parallel and perpendicular to the layers 
z = const), satisfies £oo,||/£oo,j_ = (b y /b z )~ 1 ^ 2 . The rescal- 
ing z' = (b z /b v )~ 1 / 2 z maps the large-scale continuum field 
theory of this film on a primed system with B oc 1 and film 
thickness L' = (b z /b y )^ 1 ^ 2 L. Writing the L-dependent part 
of the finite-size free energy per cross-sectional area L d_1 
and ksT at the bulk critical point as / L 4 " 1 , we can in- 
troduce a Casimir amplitude for the weakly anisotropic 
system, where p indicates which (scale-invariant) boundary 
conditions hold on sufficiently large scales. This could be 
anyone of those considered in Ref. 49, namely periodic, an- 
tiperiodic, Dirichlet-Dirichlet, Dirichlet-special, and special- 
special), as well as ++, , and H — boundary conditions. 

Upon expressing A^ p '/L d_1 in terms of L', we see that A^ 
is related to the Casimir amplitude A^ p > of the transformed 
(isotropic) system via 



= (J t |/J±r (d - 1)/2 A^ . (6.7) 



This relation was obtained for the special case of antiperi- 
odic boundary conditions (AP) in a recent paper by Dantchev 
and Griineberg 103 who determined A^ P ' in the large-n limit 
by the exact solution of a mean spherical model. Our rea- 
soning shows that it holds more generally and follows from 
simple considerations. 
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APPENDIX A: REPRESENTATIONS AND CONVENTIONS 

We define Fourier transforms with respect to time and the 
position vector's component y along the film by 

ct>(x,t) = ct>(y,z,t)= f c/>(y,z,u)e- iu ' t 

= f ' 4>faz,t)e*™ = f </>(p, z, w) e^-f-"*) ,(A1) 

where we employ the short-hand notations 

duj f f d d ~ x V 



2tt 



i-i (2tt) 



d-l 



(A2) 



Note that we do not introduce separate symbols for a func- 
tion such as <p(x,t) and its Fourier transforms <j>{p,z,t), 
<f>(p,z,u>), and <p(p,z,oj); which quantity is meant should 
be clear from the arguments of these quantities. 

When defining Fourier transforms of multi-point response 
and cumulant functions that are invariant under translations 



Hi — > Hi + Ho parallel to the boundary planes and time trans- 
lations ti — ► t i + to, we mean by the respective Fourier trans- 
forms the coefficients of the momentum and frequency con- 
serving factors (27r) d_1 S(^2 i Pi) and 27r<5(^. u>i), respec- 
tively. For example, the Fourier transforms R(x; x'; oj) and 
R(p; z, z';u>) of the free response propagator R(x, t; x' , t') 
satisfy the relations 

R(x,t;x',t') = [ R(x:x';Lu)e^ lul{t - t " ) 

J ijj 

R(p; z, z'; u) e^*"*') e ip < v ~ v '> . (A3) 
APPENDIX B: FREE RESPONSE PROPAGATOR 

To determine the free response propagator Rl(p', z, z'\ ui) 
for general non-negative values of c±, C2, and f , we must solve 
the differential equation 

[-iuj + A(f + p 2 - d 2 z )]R L (p; z, z'\ u) =6(z-z') (Bl) 

with the boundary conditions (3.7). Two linearly indepen- 
dent solutions of this Sturm-Liouville differential equation are 
e ±Ku,z^ w }j ere fc^ was defined in Eq. (3.10). From them we 
can construct the two linear combinations 



Ui(z) = k u cosh(K^z) + b\ sinh(K w z) 



(B2) 



and 



U 2 {z) = cosh[K w (L - z)] + c 2 sinh[K w (L - z)} (B3) 

that fulfill the boundary conditions on the boundary planes 
z = and z = L, respectively. The Green's func- 
tion \R L (p;z,z';uj) is given by -Ui(z < )U 2 (z > )/Wi 2 , 
where z < = min(z, z') and z > = max(z,z'). 65 The nor- 
malization constant W\ 2 is fixed by the jump condition 

\[d z R L YZi-l = - 1 - This y ields the Wronskian 



EM*') U 2 (z>) 
U{(z') U 2 (z>) 



(B4) 



which in our case is independent of z' . A straightforward cal- 
culation then gives 



\R L (p;z,z';u\ci,c 2 ) 

q., s [ci sinh^z) + k^ cos1i(k w z)] [c 2 sinh[«; w (L — z')] + k u cosh[« w (L — z')] 

«2>(ci + c 2 ) cosh(K w L) + k^k^ + c\c 2 ) sinh(K w L) 
+ (z <-* z'), (B5) 

where we have explicitly indicated the enhancement variables Cj on the left-hand side for clarity. In the special case ci = c 2 = 0, 
this simplifies to 

\ z? I < in n\ of i \ cosh (^) cosh[^(L - z')\ 

A Rl{P\ z, z ; w|0, 0) = 9{z - z) : . h (z <-»■ z ) . (B6) 

smiil Lj J 
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The free response propagator for periodic boundary conditions can be determined in a straightforward fashion by performing 
the summation in Eq. (3.11) using Eq. (3.10). This gives 

2k u smh(K w L/2) 



APPENDIX C: CALCULATION OF FEYNMAN GRAPHS 



1. Layer susceptibility for periodic boundary conditions 



The first one-loop graph of \zz shown in Eq. (5.29) in- 
volves the static propagator G^ e ^(a;|f) = °X R 1 ^^ {x , uj = 
0|f ) at zero separation x = 0. To compute it at f = 0, we add 
and subtract the zero-mode contribution. Since this vanishes 
(in dimensional regularization) when f = 0, we may simply 
replace G^"^(a;|0) by G^ el \x\0). It is then convenient to 
use the latter's representation (3. 1 1) in terms of image contri- 
butions, with G^ } (a:|0) = 2~ 2 7r- d / 2 T(d/2 - l)x 2 - d . The 
term with m = vanishes since G^ } (O|0) = in dimen- 
sional regularization. The summation of the remaining terms 
is straightforward, giving 



GT ) (« = 0|f = 0) = ^_g c(d _ 2) 



(CI) 



Since the k = component of this quantity vanishes when 
f = 0, we can directly use this result to compute the one- 
loop graph with the dashed blue line of Eq. (5.29). Upon 
performing the required z-integration of the external legs 



R£ (Pj z i w |0) 2 > one obtains 



* \ 
I I 
A O < > 9 < < O 

z z 



n + 2 u V{d/2 - 1) ((d - 2) 
3~ 2 2ir d / 2 L d - 5 
h^L + sinh(LK w ) 



8k^L 3 sinri (k wj L/2) 



(C2) 



The remaining graph of Eq. (5.29) involves the dotted red 
line, which is L^ 1 G^ X '(y = 0|<5f| pel ' ) ) in position space. 
Hence we have 



a o « < O 

z z 



kajL + sinh(K w L) 



8^3 L 3 sinh 2 (K w L/2) 
where 5?^' is the 0(u) shift given in Eq. (4.1). 



(C3) 



2. Surface susceptibilities for sp-sp boundary conditions 



The analogs of the graph (C2) for Xu ^ an d 
Xi£ involve the static propagator G^~^ \x, x'\0) = 
A ( x i x 'i w = 0|0) at equal positions x = x'. Since 



its zero-mode contribution vanishes at f = 0, it agrees with 
G^' N \x,x\0). We use its representation (3.8) and take 
into account that the m = contribution of the first sum 
vanishes at f = 0. The remaining terms can be summed in a 
straightforward fashion to obtain 



Gf N) (*,*|0) - T[d,2 - l) 



[2 CO* -2) 



2 d 7r d/2 L d-2 

+ C(d - 2, z/L) + C(d - 2, 1 - z/L)] . 

(C4) 

Note that the two generalized Hurwitz zeta functions 104 
C(2 - e, s) and C(2 - e, 1 - s) behave as s £ " 2 and (1 - s) e ~ 2 
at small values of their arguments s = z/L and 1 — s, respec- 
tively. They contain the uv singular contributions — e _1 #'(s) 
and -e- l 5'{\ - s). 

To show this more clearly and to determine the Laurent ex- 
pansion of the graph in question to order e°, we proceed as 
follows. We transform to the variable s. The z-independent 
part of Eq. (C4) leads to contributions of the graphs that can 
be expressed in term of the integrals 



Ij,k{ku,L)= I ds /j(s;k w L) f k (s;k u L) , j,k = l,2, 
Jo 

(C5) 

with 



h{s;x) 



cosh[(l — s)x] 
x sinh x 



and 



/ 2 (s;x) = /x(l - s;x) 



(C6) 



(C7) 



where Lf^z/L; h w L)/\ = i?f~ N) (P; *, 0' " 1)^5 w|0) rep- 
resents an external leg attached to 25.,. 

The required integrations are elementary, giving 

h,i{x) = h,2{x) = , Q , , 9 (C8) 



and 



h,2(x) = h,i(x) 



4x 3 sinh 2 x 



1 + x coth x 
2x 3 sinh x 



(C9) 



The contribution produced by the term oc £(2 — e, s) in 
Eq. (C4) involves the integrals 



ds£(2-e,s) /,-(>; x)/ fc (>;x) (C10) 



with (j, k) = (1,1) and (1,2). The analogous contribution 
implied by the term oc £(2 — e, 1 — s) can also be expressed 
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in terms of these integrals, as can be seen by making a change 
of variables 1 — s — > s and using Eq. (C7). 

To compute the Laurent expansion of these integrals, we 
substitute 

C(2-e,s) =s e ~ 2 + C(2-e,s + l) (Cll) 

and decompose each one of them into a sum of the respective 
two integrals 

POO 

S jik (e; x) = / dss e - 2 x) f k (s; x) 9(1 - s) (C12) 



and 



i?,- fe (e;x) = / dsC(2-e,«+l)/j(a;^)/fc(*;x). (C13) 
Jo 

The latter integrals are regular at e = 0, and hence can 
be expanded as Rj^(e\x) = i?j ; fc(0;x) + O(e). The in- 
tegrals 5j,fc(e; x) may be viewed as the results of applying 
the distribution denoted s^ 2 in Ref. 100 to the test functions 
h(s) = fj(s; x) /fc(s; x) 9(1 — s). The Laurent expansion of 
this distribution is well known. 100 It reads 

s f + - 2 = ^-5'(s) + S - 2 + 0(e), (C14) 

where the generalized function s+ 2 is defined by (cf. the Ap- 
pendix of Ref. 7) 



(0), , _ IE ~ Chi(2x) + ln(2x) - 1 

^12 \ K ) — 

x smn x 

xShi(2x) - 1 
x 2 sinh(x) tanh x ' 



(C20) 



and 



^o )(x)= xShi(2x)-coah x ; 



x 2 sinh 2 x 

where Chi and Shi are the hyperbolic cosine and sine integrals, 
respectively. 

It is useful to introduce the combinations 

S(x) = S<°>(x)+S<°>(x) 



2 coth(/i) 



[j E - 1 - CM(2k) + 1ii(2k) 



coth(K)Shi(2K)] 



1 + 2csch^x 



, (C22) 



csch 2 x ' 



R(x) = i?i,i(0;x)+ #2,2(0; x) 

ds (s + 1){2 + cosh(2sx) 

Jo 

(C23) 



2^ 2 Jo 
+ cosh[2(l-s)x]}. 



dss+ 2 h(s) = / dss" 2 [/i(s) - h(0) - s/i'(0)] 
) Jo 

/oo 
rfss -2 [/i(s) - h(0)} . (C15) 

Utilizing these results, one arrives at 



S,- fc (e; x) = 5,) fe 1} (x) e" 1 + S^(x) + O(e) (C16) 



(0), 



with 



= /j(0;x)/ fe (0;x) + /,(0;x)/a0;x) 



1 



2 coth(x) , j = k = 1 , 

--^csch(x), j^fc, 
x 

0, j = * = 2, 



(C17) 



and 



5 $W = {/ J ( S ;^/ fc ( S ;^-/ J (0;x)/ fc (0;x) 
Jo s 

- « [/j(0; x) / fc (0; x) + /,(0; x) /£(0; x)]} 
-/ 3 -(0;x)/jfc(0;x), (C18) 

where the prime on f'j(s; x) means a derivative with respect 
to s. 

The latter integrals can be performed in a straightforward 
fashion using MATHEMATICA. 105 This yields 



Su(x) = [IE - Chi(2x) + ln(2x) - 1] 

^ x cosh(2x) Shi(2x) - 1 
x 2 sinh 2 x 



2 coth x 



(C19) 



and the function 



i? 2 (x) = 2i? lj2 (0;x) 



2csch 2 x ' 



ds{V (1) (« + 1) 
x cosh(sx) cosh[(l - s)x]} . (C24) 

Then our results for the graphs involving the dashed blue line 
can be written as 



I I 



n + 2uAr(l- £ /2)L 1+e f 2coth(«; w i) 



3 2 (47r) 2 - e / 2 



Ku,Le 



■ h,x{ku,L) + S(k u L) + R(k u L) + O(e) \ (C25) 



and 



I I 
A (8 < W < 8 



»1 ®2 

n + 2 u\ T(l - e/2) L 1+£ ( 2 csch(K w i) 



3 2 (4tt) 2 ^/ 2 



k,.,Le 



y hA^L) + 2Si 0) 2 (k u L) + R 2 (k UJ L) + O(e) 



(C26) 



21 



The remaining graphs, involving the zero-mode propagator 
at the shifted bare temperature <5-f| sp sp ', are given by 

= n + 2 u A^i [(5f (sp-sp) ](d _ 3 ) /2 h j{kuL) (C27) 

While the latter graph is uv finite, the previous two contain 
pole terms. Recalling the 0(u) result 7,66 ' 67 Z\ = 1 + u (n + 
2)/3e + 0(u 2 ), one sees that they cancel with the contribution 
(Z^ZT 1 - l)A/i- 4 ® — <-h». 

In our calculation of the zero-momentum limits of the in- 
verse oj = susceptibilities Xn jj an d X\l w trie asymptotic 



behaviors of the functions S, R, i?2, and S 1 }^ as >f ^ 
are needed. Using the known limiting forms Chi(xr) = 
Je + In >f + 0(>*r 2 ) and Shi(xr) = x + 0(>*r 3 ), one finds 
that 

S{x) = 2S[° 2 } (x) + 0{k- 2 ) = -2>^ 4 + 0(k- 2 ) (C28) 

for small k. To determine the asymptotic forms of the func- 
tions R(k) and R2(x), we expand the integrands in the inte- 
grals of Eqs. (C23) and (C24) in powers of x. Performing the 
s-integral for the lowest-order term and expanding the prefac- 
tors leads to 

R{x) = R 2 {x) + 0(x- 2 ) = 2k- a + 0{x- 2 ) . (C29) 
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